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Abstract 


The Weyl double copy is a procedure for relating exact solutions in biadjoint scalar, gauge 
and gravity theories, and relates fields in spacetime directly. Where this procedure comes from, 
and how general it is, have until recently remained mysterious. In this paper, we show how 
the current form and scope of the Weyl double copy can be derived from a certain procedure 
in twistor space. The new formalism shows that the Weyl double copy is more general than 
previously thought, applying in particular to gravity solutions with arbitrary Petrov types. We 
comment on how to obtain anti-self-dual as well as self-dual fields, and clarify some conceptual 
issues in the twistor approach. 


1 Introduction 


The double copy is by now a highly-studied procedure for turning solutions in a (non-)abelian 
gauge theory into gravitational counterparts. It originally arose in the study of perturbative scat- 
tering amplitudes m], where it had a string theoretic motivation at tree-level. Although the 
double copy remains conjectural, there is by now a huge amount of evidence for its being exact in 
a variety of (supersymmetric) theories, including at the quantum level, and to all orders in per- 
turbation theory in certain kinematic limits (aHa). A related body of work has extended the 
double copy to classical solutions. The first work to appear was the Kerr-Schild double copy of 
ref. [44], which concerned certain exact (albeit algebraically special) solutions of the Einstein equa- 
tions, and demonstrated the existence of well-defined counterparts in gauge and biadjoint scalar 
theory. A second exact procedure is the Weyl double copy of ref. [45], which is more general than 
the Kerr-Schild approach, although equivalent to the latter where they overlap. However, it is 
still apparently restricted to algebraically special solutions, albeit ones that may be interesting for 
astrophysical purposes (see e.g. refs. for follow-up studies). A key feature of exact classical 
double copies is that the equations of motion in each theory turn out to be linearised, with no 
non-linear corrections. Other classical double copy techniques have been developed that can in 
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principle go beyond this (82}/119], but at the price of proceeding order-by-order in perturbation 
theory, such that one loses an exact understanding. One might thus hope that further scrutiny of 
the Kerr-Schild and Weyl double copies - including ascertaining their limitations and scope - will 
yield an underlying explanation of where the double copy comes from, and an understanding of 
how to apply it to arbitrary solutions. Given the fact that the Weyl double copy is more general 
than the Kerr-Schild double copy, it is sufficient to focus on the former. 


Recently, a derivation of the Weyl double copy has been given [120], using ideas from twistor 
theory [1211123] (see e.g. refs. for pedagogical reviews). This explained why the Weyl 
procedure has the form it has, as well as suggesting that it is more general than previously thought. 
In particular, the original Weyl double copy of ref. (see also ref. [128}) applies to gravitational 
solutions of Petrov types D and N onlyf} However, ref. found an example outside these classes, 
and also argued that the Weyl double copy should extend to general conformally flat spacetimes, 
thus formalising preliminary remarks in this regard (see also for a recent discussion of 
the classical double copy in curved space). The aim of the present paper is then twofold. Firstly, we 
will fill in the details of the brief ref. [120], providing full details of how to carry out the appropriate 
calculations. Secondly, we will extend previous results to anti-self-dual as well as self-dual fields, 
as well as to arbitrary Petrov types. 


The structure of our paper is as follows. In section |2| we will introduce salient details regarding 
spinors and twistors, as well as reviewing the Weyl double copy of ref. [45}. In section [3] we will 
present our twistor-space formalism for obtaining spacetime Weyl double copy formulae, going 
beyond the preliminary results of ref. [120]. We summarise our results and conclude in section [4] 


2 From spinors to twistors 


The Weyl double copy relies on the spinorial formalism of General Relativity and related theories. 
Although this is textbook material (see e.g. ref. in addition to the above references), this 
formalism is not necessarily known to all researchers working on the double copy or beyond, and 
the same can certainly be said about twistor theory. We will thus review key concepts in this 
section, in order to make our presentation self-contained, and also to set up crucial notation needed 
for the rest of the paper. 


2.1 The spinorial formalism 


Our first introduction to GR typically uses the language of tensors and four-vectors. However, 
an alternative formulation exists, in which all equations are expressed in terms of two-component 
spinors m^ = (n, nt), and their higher-rank generalisations. Spinor indices may be raised and 
lowered according to 


TA =EABn?, nP =r4e^P, (1) 


where €4p and <4" are the two-dimensional Levi-Civita symbols defined such that 


cape’? = 65, en =1, (2) 


4We review the Petrov classification below, after introducing the appropriate language. 
°Our conventions follow those of refs. 125), where we have chosen an appropriate spin basis to define €o1. 


where oa denotes the Kronecker symbol. Given a spinor 74, one may also consider its complex 


conjugate r^, where primed indices may be raised and lowered analogously to eq. (i), but with 
the symbols ^B" and eyg’, such that 


si / 1 
expel? =G, egy =1. (3) 


These operations extend to objects with any number of spinor indices. Furthermore, there is a 
remarkable simplification of the structure of these higher-rank quantities, which ultimately follows 
from the fact that there are only two values that can be held by each spinor index A or A’. Let us 
introduce the notation 


1 
P(A142...An) = n! 5 PA (1) Ac(2)--Ac(ny? (4) 


where the sum is over all permutations o of the index labels (1,...,). That is, eq. constitutes 
the fully symmetric combination of spinor components of an arbitrary rank n spinor (with suitable 
generalisation to primed indices). Then any multi-rank spinor can be decomposed into a sum of 
terms, each of which involves symmetric spinors, multiplying Levi-Civita symbols. We will see 
explicit examples of this shortly. Another nice property is that any symmetric spinor factorises 
into a symmetrised product of spinors e.g. 


S4p..c =Scap..c) = SAB..c = APB--- Yc): (5) 


The individual spinors {a,4, 3g,...} are associated with null vectors in spacetime, referred to as 
principal null directions of Sap..c. To see this, one may note that any tensorial quantity can be 
translated into the spinorial language using the so-called Infeld-van der Waerden symbols {oh AD 
which may be chosen in Cartesian coordinates as follows: 


1 1 0 i 1 0 1 / 
0 _ „AA 1 = _ „AA 
oy = ( 0 1 ) = 00 ; TAA = ( 1 0 ) = 0i , 
I 0 =i , il 1 0 , 
2 = —_ „AA 3 = _ „AA 
TAA = V2 ( i 0 ) = —02 , GAA V2 ( 0 —1 ) a3 3 (6) 


For a 4-vector this gives 


a 1 f/f Wt YU-ive 
veiw =e tak V — V3 , (7) 
where the determinant of the matrix on the right-hand side is 
1 
det (Vaoaa) = 5 (Vo)? — (Vi)? — (Va)? — (V3)?) . (8) 


We may recognise this as being proportional to the norm of the 4-vector, such that the determinant 
vanishes if Va is null. By standard linear algebra arguments, this then implies that the matrix must 
factorise i.e. 

VaV% =0 =} Vi644 = TATA, (9) 


where m4’ = (7,4)* given that the matrix in eq. is clearly Hermitian. Conversely, given any 
spinor 74, we may construct a matrix MAy = TATA, which in turn corresponds to a null 4-vector 
in spacetime. In particular, each of the so-called principal spinors appearing in the decomposition 


of a general symmetric tensor (eq. (5)) can be associated with a principal null direction in spacetime. 


A given solution to Einstein’s field equations in GR will have a corresponding Riemann tensor 
Ragys- One may translate this into the spinor language as above, and then decompose it into various 
symmetrised spinor parts, where some further simplifications arise due to the known symmetries 
of the Riemann tensor itself. The result turns out to be 


Rasy > Raa BB'CE' DD = VaBcpea’B'€c'p! + Val BIC’ D!EABECD 
+ BABE D'EA B'ECD + ®A'BICDEABEC'D! 


+ 2A(€aceBDEA'B'EC'D! + EABECDEA'D'EBIC’); (10) 


where all spinors appearing on the right-hand side are fully symmetric. The spinors ®4g,4/p and 
®4/p AB are directly related to the trace-reversed Ricci tensor Rag, and A is directly proportional 
to the Ricci scalar R (see e.g. ref. (126). We will be concerned with vacuum spacetimes, so that 
these quantities all vanish by the Einstein equations. We are then left with the only contribution to 
the curvature that is present in free space, which in the usual tensorial formulation of GR is known 
as the Weyl tensor, and denoted Cag 5. From eq. (10), we thus have the spinorial identification 


Capys + YABCDEA B'E D + Vai BCD! EABECD- (11) 


If we are working in Lorentzian signature such that the spacetime is real, Y y p/crp must simply 
be the complex conjugate of V4gcp. More generally, Vagcp and Yy p'ap are the anti-self-dual 
and self-dual parts of the Weyl tensor respectively. That is, they are respectively projected out by 
the conditions 

"Capys = FiCabys, (12) 


where the dual Weyl tensor is defined by 


1 
*Capyé = z labor On yâ: (13) 
The dynamics of the Weyl tensor is constrained by the Bianchi identity for the Riemann tensor, 
which can be shown to lead to the following conditions: 


V44'Uapcp =0, V“ Uy pcrp = 0, (14) 


where V44" = vo is the appropriate spinorial translation of the spacetime covariant derivative. 
Given its role as part of the Weyl tensor, the spinor V 4gcp is usually referred to as the Weyl spinor. 


Let us now turn to electromagnetism, whose equation of motion involves the field strength tensor 
Fv. Using similar methods to those above, one may show that the spinorial decomposition of the 
field strength is as follows: 


Fug > FAA'BB' = ỌABEA'B' + ỌA'B'EAB, (15) 


where the symmetric spinors ¢4g and ġy'p are the anti-self-dual and self-dual parts, respectively 
projected out by 


* a * 1 OT 
Fag = FiF op, Fag = gopor : (16) 


The Maxwell equations then imply 
VAi bag = 0, Vo gyp = 0. (17) 


We may note that eqs. and are both special cases of the general spinorial equations (see 


e.g. ref. [124]) eee 
V4“ bap..c=0, VA" yp. ao =0 (18) 


where ¢48...¢ is assumed symmetric, with n indices. These are known as the massless free field 
equations, as they are indeed associated with massless and non-interacting fields in a vacuum 
spacetime. The spin of the field is given by the number of spinor indices divided by two, which 
matches the two cases described above: eqs. and contain spinors with four and two spinor 
indices, describing the spin-2 graviton and spin-1 photon respectively. If we restrict to solutions 
of positive frequency, the spinors ¢4/p’_.cr and d,4g..c represent states of positive and negative 
helicity +n/2 respectively (in units of ñ) (| Note that we have not yet stated which spacetime we are 
working with, which above corresponds to the fact that V4" is the covariant derivative associated 
with a potentially curved spacetime. In what follows, we will only need to consider eqs. in 
Minkowski spacetime, although we will not necessarily work in Lorentzian signature. Furthermore, 
we may also work in complexified Minkowski space, whose line element is 


d3? = dt? — dz? — dy? — dz”, t,2,y,z€C (19) 
n.b. the complexified coordinates, but not their complex conjugates, appear. 


An immediate use of the spinorial language is that it allows us to classify different types of solution 
in electromagnetism and gravity. Above, we noted that any symmetric spinor can be factorised 
into 1-index principal spinors. These may be degenerate, and the various different patterns of 
degeneracy allow for a classification of different solutions. For the electromagnetic field strength 
spinor, one has: 

PAB = APB), (20) 


and there are then two different “types” of field strength spinor: (i) those with distinct null direc- 
tions (aa & a); (ii) those with a degenerate null direction, so that a4 « 84. The latter give rise 
to null electromagnetic fields, given that the field strength spinor then satisfies 


banger? = A PP pangon = Ce Faaapacap = 0. (21) 
For the Weyl tensor there are more possibilities. In general one has 


YABCD = AbBYCOp); (22) 


and one can classify the different types of Weyl spinor by their pattern of degenerate null directions. 
If they are all different, this is called a spinor of type {1,1,1,1}. Gradually making more of the 
null directions degenerate leads to types {2,1,1}, {3,1} and {4} (n.b the symmetrisation of the 
l-index spinors in eq. means that the order of the principal spinors does not matter - only 
their degeneracy). A fifth possibility is that there are two distinct pairs of null directions, denoted 
as {2,2}. Finally, the Weyl spinor may be zero, which is written as {—}. This reproduces the 


°In optics parlance, positive and negative helicity correspond to right-handed and left-handed circular polarisations. 


Weyl type | Petrov label 
{1,1,1,1} 
12.117 
{3,1} 
{4} 
{2, 2} 
{=} 


Table 1: Different types of Weyl spinor classified by: (i) the pattern of degenerate principal null 
directions; (ii) the equivalent Petrov type. 
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well-known Petrov classification of GR. solutions, which was first derived using tensorial methods. 
That formalism used different labels to those used here, and we summarise these in table|1} We 
will use the Petrov labels from now on. 


Having summarised spinorial notation and methods, we can now state the Weyl double copy dis- 
cussed in the introduction, which was first presented in ref. [45]. Given an electromagnetic field 
strength spinor @4pB, one may construct a Weyl spinor according to the rule 


1 
Wascp = gP(ABPCD)) (23) 


where S is ascalar function. This procedure was argued to hold for arbitrary type D vacuum space- 
times in ref. , where the scalar S could then be found in particular examples by matching both 
sides of eq. All of these solutions have the property that they linearise the Einstein equations. 
Thus, the corresponding field strength and Weyl spinors may be taken to satisfy eqs. 17) in 
Minkowski space i.e. such that V4" corresponds to a flat-space derivative. The rule of eq. (23) is 
given only for the anti-self-dual part of the Weyl tensor. It is straightforward to write down the 
appropriate generalisation for the self-dual part, and the above linearity property then means that 
one may superpose these solutions to obtain the complete Weyl tensor. 


Equation makes intuitive sense from our previous discussion of principal null directions. Taking 
a field strength spinor of the form of eq. (20), one obtains a Weyl spinor of form 


Wascn ~ (4bBacBp), (24) 
which is clearly of Petrov type D. This is not the only way to obtain such a Weyl tensor. As already 
argued in ref. [45], one could also define a “mixed” Weyl double copy 

Vascp = Séanĝon), (25) 
involving two different electromagnetic spinors 

aB =a4aB, bcp = BoBp. (26) 


This immediately suggests that one might be able to generalise the Weyl double copy away from 
gravity solutions of type D: by taking different patterns of principal null directions of both electro- 
magnetic spinors in eq. (25), one may obtain a variety of Petrov types, as summarised in table [| 


PAB bc p | Petrov type 
aaBbB | ycop 
aaBB | AcYD 
asap | Boyp 
aaag | acp 
QAQB | ACAD 
asap | BcBp 
aABB | acBp 


Skok prea gH 


Table 2: Various combinations of principal null directions for two electromagnetic spinors lead to 
Weyl tensors of different Petrov types under the mixed Weyl double copy of eq. (25). 


We have yet to show whether such a double copy is meaningful, although we will see later that this 
is indeed the case, provided one restricts to linearised level in some cases. 


Since its inception, the Weyl double copy has been used to study certain topologically non-trivial 
electromagnetic solutions (“Hopfions” ) [130], and a tensorial formulation of the relationship between 
the electromagnetic field strength and Weyl tensor has been presented in refs. (69\[71), allowing a 
generalisation to higher dimensions. An alternative approach to the self-dual double copy which 
also has links to spinorial ideas may be found in ref. [64]. Recently, a study of the Weyl double 
copy properties of type N solutions has been presented in ref. [128]. 


2.2 Twistors 


In what follows, we will be concerned with deriving the Weyl double copy using ideas from twistor 
theory (121}/123}. This is a subject with an illustrious history spanning more than half a century, 
and we cannot do justice to all of its many developments here, but will concentrate on those aspects 
which are crucial for what follows. Excellent reviews may be found in refs. (125}127]. We may 
start by defining twistor space T as the set of solutions of the twistor equation 


Vn) =0, (27) 
whose general solution in Minkowksi space is 
04 = w4 — iA’ ny. (28) 


One may thus associate solutions of eq. with four-component objects (“twistors”) containing 
a pair of 2-spinors: 
Le = (w^, ra) , (29) 


More explicitly, the index notation on the left is such that 
(Z0, Zt, Z?, Z3) = (w, wt, ny, my). (30) 


The “location” of a twistor in Minkowski space is defined to be the region in which its associated 
spinor field Q4 vanishes. From eq. (28), this implies the incidence relation 


wit = A’ ny, (31) 


A 


where 74’ is the spinorial translation of a point in space-time. Note that this condition is invariant 


under simultaneous rescalings 
A A 
we > hw, mwa Ata, AEC, (32) 


so that twistors obeying the incidence relation are defined only up to an arbitrary complex scale fac- 
tor. They thus correspond to points in projective twistor space, or PT. If we consider a fixed point 
x^^ in eq. (31), this defines a line in PT. To see why, note that the constraint of eq. reduces the 
four complex parameters in a general twistor to two, and a further complex parameter is removed 
by going to the projective space. We therefore see that the relationship between Minkowski space 
and PT is non-local|"} A complex line (with a point at infinity) can be mapped to a Riemann sphere. 


By considering the conjugate twistor equation 


VEAP) =0 (33) 
whose general solution is 
AA = u” big d, (34) 
one may define dual twistors 
Wa = (a " a (35) 


Projective dual twistor space is denoted by PT*. Similarly to twistors, the location of a dual twistor 
in spacetime is defined by the region where the associated spinor field A4 vanishes, leading to an 
incidence relation 

u” = ied" ja. (36) 


There is an inner product between twistors and dual twistors, defined by 
Z“Wa = wir, T un TA, (37) 


where we have defined components as in eqs. B5). These various objects have a number of very 
nice mathematical properties. For example, one may show that the generators of the conformal 
group act linearly on (dual) twistors, and that the inner product of eq. is conformally invari- 
ant. This explains the ubiquity of twistor techniques in modern research on scattering ampitudes, 
given that there are many situations in which (dual-) (super-) conformal invariance can be made 
manifest, particularly in more symmetric theories such as M = 4 SYM. 


Another useful result from twistor theory is that there is a known correspondence between solutions 
of the massless free field equations of eq. in Minkowski spacetime, and certain contour inte- 
grals in projective twistor space (122). This is called the Penrose transform, and involves twistor 
functions [Ë| f(Z°) (i.e. not involving the conjugates 7“). Given such a function with Z® defined 
as in eq. (29), we may define the integral 


1 / 
pap. alt) = Iri $ ngdr” TATB .. TC | Prf (Z), (38) 


TLikewise, points in PT represent two (complex) parameter surfaces called a-planes in complexified Minkowski 
space, or a null ray through a given point in real Minkowski space. 
8The functions f(Z%) are actually representatives of cohomology classes, as we describe below. 


where the symbol p, denotes that we must restrict to the line in PT corresponding to the spacetime 
point x^^. The contour I for this integral is defined on the related Riemann sphere, and is well- 
defined only if the various poles appearing in f(Z°) can be separated from each other. We will see 
explicit examples of functions and contours later on. Note that, for eq. to make sense as an 
integral in projective twistor space, the integrand (including the measure) must be homogeneous 
of degree zero under rescalings ty + Ary (or Z% —+ AZ). This in turn implies that the function 
f(Z*) must have degree (—n — 2), where n is the number of indices appearing on the left-hand side. 
Clearly the field on the left-hand side is symmetric in { A’, B’,...,C’}, and thus could potentially 
represent a solution of eq. (18). The proof that it does indeed do so is relatively straightforward. 
First, one notes 


Vonet Z] = yli nany) 
aoe E | (39) 
Then from eq. we have 
VDD B... = TU -< -TOT D! Px El l (40) 
The right-hand side is symmetric in A’ and D’. Thus, contracting both sides with ^P” yields 
Vi papcr =0 (41) 


as required. Note that there is considerable freedom on the right-hand side of eq. (38), due to 
being able to move the contour, and also being able to modify the function f(Z°) with additional 
terms that vanish upon performing the contour integration. To examine this in more detail, note 
that any contour will separate the Riemann sphere into two parts, such that any singularities are 
associated with two regions on either side of the contour. Without loss of generality, we denote 
this as in figure [1] with the contour at the equator of the sphere, and the singularities confined to 
regions N and S in the northern and southern hemispheres respectively. We must get the same 
answer for the integral of eq. if we take the contour to enclose wither N or S (up to a change 
of sign for differing orientation of the contour). However, we are free to modify the twistor function 
according to the equivalence relation 


F(Z°) > F(Z°) + fn (Z°) + fs(Z%), (42) 


where fy(Z°) (fs(Z)) contains singularities only in the northern (southern) hemisphere. In carry- 
ing out the contour integral for either of fv,9(Z%), one may simply choose to enclose the hemisphere 
that is free of singularities, thus obtaining zero. 


The above remarks make clear that the Penrose transform map from twistor space to spacetime 
is many-to-one. In more mathematical terms, this may be described in terms of sheaf cohomology 


©The fact that only the derivative with respect to w? arises in eq. (39) and not its conjugate, is a consequence of 


f(Z°) being a representative of a cohomology class. See ref. for a discussion of this issue from a different point 
of view. 


Figure 1: A contour I divides the Riemann sphere into two hemispheres, with singularities in two 
regions N and S. 


groups (see e.g. refs. for a review), where the notation H! (PT, O(m)) represents, roughly 
speaking, the equivalence class of holomorphic functions of homogeneity m in projective twistor 
space (avoiding the regions N and S), modulo the redefinitions of eq. (42). The above Penrose 
transform then amounts to an isomorphism between massless fields of helicity n in spacetime, and 
the cohomology group me (PT, O(—n — 2)). 


Above, we have given the Penrose transform for primed fields. A similar formula may be written 
for unprimed fields, if one instead uses holomorphic functions defined on projective dual twistor 
space: 


ie = = $ E E A A (43) 


Analogously to eq. (38), this constitutes an isomorphism between massless fields in spacetime with 
helicity —n, and the cohomology group H!(PT*,O(—n — 2)). In early works on twistor theory, 
it was seen as unsatisfactory that both twistors and dual twistors were needed to define fields of 
arbitrary helicity. This in turn led to the alternative procedure for negative helicity fields 


— 1 E' 0 Q 
WaAB...c(2) => On f rwan Px so Bie tee A ot 4 ) , (44) 


in which the spinors w^ are defined in eq. (29). In this formula, the twistor function must now 
have homogeneity (n — 2) for consistency. In terms of the sheaf cohomology groups mentioned 
above, one has an isomorphism between helicity —n massless fields in spacetime and elements of 
the cohomology group H!(PT, O(n — 2)). From now on, we will loosely use the term function to 
imply a given representative of the appropriate cohomology class, returning to a discussion of the 
latter in section 


Each of eqs. is defined for Minkowski space only. In gauge theory or gravity, they thus 
give rise to field configurations that satisfy the linearised Yang-Mills or Einstein equations. This 
does not bother us for the Weyl double copy, all previous examples of which have been linear, but 
exact, solutions. However, we must bear this in mind when considering generalisations. Full non- 
linear generalisations of the Penrose transform exist for the (anti-)self-dual sectors of Yang-Mills 


Tn practice, one takes this isomorphism on a suitable open subset of PT e.g. that corresponding to positive or 
negative frequency fields in spacetime. 
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theory and gravity [133]. Also, even the linear Penrose transform above is more general than it 
seems. Conformal invariance is manifest in twistor space, so that we expect that the fields obtained 
from eqs. and (44) can also be transformed to an arbitrary conformally flat spacetime. In fact, 
the massless free field equation is known to be conformally invariant [124], so this is indeed the 
case. There is a subtlety, however, regarding scalar fields. From eq. (38), we expect these to be 
given by contour integrals of the form 


277 


b(n) = — $ redr” [psf (Z°). (45) 


We can longer form an equation like eq. (18), as there is nothing for either of the spinor indices in 
the covariant derivative V4“ to contract with. However, one can instead show that eq. is (in 
conformally flat spacetimes) a solution of the conformally invariant wave equation 


( +E) o=o, (46) 


where R is the Ricci scalar. 


The above-mentioned many-to-one property of the Penrose transform of eq. means that it is 
not possible to write a unique inverse transform that fixes a twistor function f(Z®) from a given 
spacetime field. However, there are some tricks for formulating representative twistor functions for 
spacetime fields possessing certain properties. A particularly useful one is the observation made 
in ref. (125), relating the poles of twistor functions to principal spinors in spacetime. First, note 
that the factorisation property of symmetric spinors means that if a given n-index spinor has a 
k-fold principal spinor £x, it will vanish if contracted with (n — k + 1) factors of £x, but not if 
only (n — k) factors are contracted. To see this, we may write such a spinor’s decomposition into 
principal spinors as 


QABE! = AEB -EC ap Be... YF); (47) 
ee n m 
k factors (n—k) factors 


so that contracting this with (n — k) factors of £^ gives 


EP'e.. E barr = [a€l [BE]... [yél&caréw---€0n, Au] = Aare’. (48) 


The various prefactors on the right-hand side are manifestly non-zero, given that none of the spinors 
{a’,, Bp... yor} can be proportional to £x, if they are distinct principal spinors. Contracting with 
a further factor of £^ immediately gives a factor of [E£] = 0, thus proving the assertion made above. 
Now consider the spacetime field defined by eq. (38), in the case where the twistor function f(Z°) 
has a pole of order m (where m < ny which occurs when mty x na. Contracting eq. with 
m factors of n” gives 


1 


në nP B n” PA'B'...C'D!...F'(L) = oa frea” [nn nre ... nr lpr f (Z. (49) 
SS a a 


m factors n indices 


“Tf the twistor function has a pole of higher order than n, enclosed by I’, then the associated spinor is not even a 
principal spinor. 


11 


The factor [77] in the integrand will kill the m**-order pole as 74: > na’, such that the contour 
integral is zero. By the above remarks, and given that this will not occur if one contracts with only 
(m—1) factors of n^’, we thus find that in the Penrose transform of eq. (38), the field y’ p'...p' has 
at least a (n — m + 1)-fold principal spinor ņa, if the twistor function f(Z®) has a single m**-order 
pole as 74’ —> na, enclosed by T. Furthermore, if this pole remains present for varying 744’, then 
the twistor function must have the general form [125] 


F(Z) = Om (Z° x(2°) (50) 


where 0,,(Z%), x(Z%) are homogeneous and holomorphic, 6,(Z%) is regular at the m-fold pole we 
are discussing, and y(Z®) has a simple zero. 


3 A twistorial derivation of the Weyl double copy 


In the previous section, we have reviewed various aspects of spinors and twistors, culminating in 
the Penrose transform of eq. (38), and its related results of eqs. 44). In this section, we show 
how these results can be used to derive the Weyl double copy of ref. [45], presented here in eq. 
(23). We will focus on the conjugate form of eq. (23), involving fields with primed indices. From 
the remarks of section [2.2] this means that we will be concerned with maps from PT to spacetime, 
rather than PT*. Changing notation for later convenience, the general (mixed) Weyl double copy 
may be written as 


1 
ỌA'B'C'D' = zastop (51) 


Here ¢ satisfies the scalar wave equation, and oy?) are two electromagnetic spinors. Then ¢,4/pB'c" D’ 
satisfies the spin-2 case of the massless free field equation of eq. (18), and thus represents a self-dual 
Weyl spinor. As discussed in the introduction, the Weyl double copy is related to previous exact 
classical double copies that have appeared in the literature where overlap exists e.g. the Kerr- 
Schild double copy of ref. [44]. It has been argued to hold for general type D vacuum spacetimes 
in ref. |45|, and also for type N vacuum spacetimes in ref. (128]. However, a number of questions 
remain: 


1. The form of eq. involves a certain product of spinorial quantities, with symmetrisation 
over the indices. Is there a deeper explanation of where this form comes from? 


2. Similarly, why should there be an exact double copy relating fields in position space? The 
original BCJ double copy for amplitudes involves products in momentum space, which 
would be expected to lead to a convolution in position space. Indeed, there are alternative 


double copy formalisms that have exactly this property 138|. 


3. Does the double copy apply to Petrov types other than type D or N? As already discussed in 
section [2.1] the form of eqs. would seem to apply for arbitrary Petrov types, although 
it is not of course guaranteed that the resulting spinorial quantity on the left-hand side will 
satisfy the spin-2 massless free field equation. 


4. Can the Weyl double copy be extended to curved spacetime backgrounds? Preliminary work 
for the Kerr-Schild double copy implies that exact results are possible in conformally flat 


spacetimes ; 
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Reference recently presented (in a very brief form) a twistor space procedure for deriving the 
type D Weyl double copy, and also gave an example of a more general solution (of Petrov type III) 
that could be obtained from eq. (51). We discuss this, with full details that are missing in ref. (120), 
in the following sections. 


3.1 Twistor space picture 


Consider a two holomorphic twistor functions FS?) (Ze) of homogeneity —4, and a further holo- 
morphic twistor function f(Z°%) of homogeneity —2. By the Penrose transform described here in 
section 2.2} these will necessarily respectively correspond to electromagnetic spinors 072) and a 
scalar field ¢ in spacetime. One may then form a product 


(Q) yay £2) 
ay _ fr(Z*) fem (2%) 
Serav.(Z°) = F(Z") (52) 


such that the function on the left-hand side necessarily has homogeneity —6, and thus potentially 
corresponds to a spacetime field solving the spin-2 massless free field equation i.e. to a self-dual 
(linearised) gravity solution [| All we have done here is use known properties of the Penrose 
transform. However, the above remarks imply that there might be some sort of relationship be- 
tween gauge, gravity and scalar fields in spacetime that corresponds to the product of eq. (52). 
By choosing suitable functions on the right-hand side of eq. (52), we can generate a particular 
spacetime relationship. We now show that for a suitable choice of twistor functions, this spacetime 
relationship is precisely the type D Weyl double copy of eq. (51). 


To find the appropriate functions, we may rely on the result of eq. (50), namely that a spacetime 
field possessing a (n — m + 1)-fold principal spinor must have a m-fold pole in twistor space. A 
type D Wey] spinor has a pair of 2-fold degenerate principal spinors, such that the twistor function 
f(Z°) must have two triple poles. This in turn implies that the function y(Z®) in eq. has two 
simple zeros, and thus corresponds to some quadratic form 


x(Z*) = QaaZ°Z", (53) 
for some constant Qag. Let us also define Om in eq. to be the simple combinatorial factor 
(54) 


for reasons that will become clear. Then, putting everything together, we may define a family of 
twistor functions 


fn(Z*) = = [Que Z2Z"] ”, (55) 


The reader may object at this point that we are multiplying together twistor functions. We return to this in 
section 

31¢ is important here that the Weyl spinor we are seeking has no more than two distinct principal spinors: more 
than two, and the contour I in the Penrose transform would be enclosing more than one pole in at least one hemisphere 
of the Riemann sphere, thus invalidating the remarks leading to eq. (50). We discuss this point in more detail in 


section 
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and our claim is that this will produce a type D Weyl tensor (for m = 3), that is related to an 
electromagnetic field strength (m = 2) and scalar field (m = 1). To show that this is indeed true, 
we may carry out the Penrose transform in each case by choosing homogeneous coordinates 


TA = (155); E EC (56) 


(n.b. we are allowed to do this, given that the integrand of the Penrose transform has homogeneity 
zero under rescalings of ty). The quadratic form in eq. is to be evaluated when the incidence 
relation of eq. is obeyed, namely for Z® = (ix^ ry, ma). Substituting eq. then implies 


x(Z°) = NH (x) (E — &)(E — £2) (57) 


in general, where €; = &(x) expresses the location of a pole in terms of the variable £, and N~!(a) 
is an overall normalisation factor. The dependence of the incidence relation on x^^ gives rise to 
the spacetime dependence of {N,€;}. The measure of the Penrose transform becomes simply 


a” dre = dé, (58) 


and we can now consider each value of m in turn. For the scalar case of m = 1, the Penrose 
transform reduces to 


1 N(x 
$= pas (2) 
ani Jr ° (E — £1) (€ — £2) 
N 1 
(x) f dé l 
2mi Jp > (E— E)E- £2) 
The contour T is defined on the Riemann sphere of £, such that it separates the poles at € = & 
and € = 9. We may thus evaluate the contour integral by enclosing only one of these poles, which 
we take to be € = é. Cauchy’s theorem then tells us that 


(59) 


dE 2ri 
= 60 
leae E&i — êz’ (60) 
so that we have 
_ N(z2) 
Next, we have the case m = 2, which gives the Penrose transform 
_ N?(x)1 (La (1, £p 
QAB = mi D! fa (E — & (E — fo) (62) 


We choose to evaluate the integral for each choice of { A’, B’} by again enclosing the pole at € = &. 
Recalling that if a function f(z) has an n‘*-order pole at z = c, the residue associated with the 
latter is 


1 i d=! 7 
res(f, c) = qy lim Fo le- oo], (63) 
one may then verify the following: 
1 en i 2; n = 0, 
an = = 4 
ini J. EEE EET” ow = 
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which in turn implies 


N*(z) N*(z) 1 N?(x) 

1 = - > TT p= in = } : 11 = — š 65 
Po'o GETAL do = vo (4 —&)3 5 (61 + &2), ba E gree (65) 

Alternatively, one may express ¢,4/p directly in terms of its principal spinors. Defining 
ax = (1,1), Bar = (1, £2), (66) 

eq. (65) is equivalent to 
TETE. ae (67) 
WP Er = &)8 AEB 


Finally, we can examine the case m = 3, which produces a Weyl spinor as follows: 


N? (x) it (1, €)ar(1, €)ar(1, €or (1, €) pe 
m a fe eee GP 


QA B'O D = 


Again using eq. (63), one finds 


6, n=V, 
3(&1 + £2), n= 4, 
1 g" = 1 2 2 = 
oni f. “EEE E T * a a = 
GEZER, n= 4. 


Substituting this into eq. (68), one finds that the result may be written as 


N? (x) 


ABCD = E = & BaP ac’ Bp!) 


N(x) | 
=|—— | ara ocd); (70) 
R =] 
where in the second line we have recognised the spin-1 massless field (i.e. an electromagnetic field 
strength tensor) of eq. (67). and used the symmetrisation property 
WV ((A4'B!)(C'D')) = VA B'O D)» (71) 
which holds for an arbitrary spinor V4/p/¢1p:.. We may recognise the prefactor in eq. as the 
inverse of the scalar field of eq. (61). Thus, our choice of twistor functions has provided a scalar, 
electromagnetic field strength and (linearised) Weyl tensor satisfying 
1 
PA'BIC!D! = gle boD): (72) 
This is precisely the (primed version of) the Weyl double copy of eq. (23). Some further comments 
are in order: 
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e Our choice of the combinatorial factor in eq. is to reproduce the normalisation of the 
Weyl double copy as given in ref. [45], and one needs such a factor whenever higher-order 
poles are present. However, the scalar function S(x) is itself only defined up to a constant 
factor in that paper, so that this is not strictly necessary. Furthermore, we are working with 
linearised field equations, so that any constant factor is possible. 


e The twistor story explains why there is a classical double copy in position space, given that the 
Penrose transform links twistor functions with spacetime fields. Furthermore, the somewhat 
mysterious form of eq. is now also explained. 


e In the original formulation of the Weyl double copy, there was no clear prescription for fixing 
the scalar function S. Here, we see that it naturally arises as the scalar field ¢ obtained from 


the m = 1 case of eq. (55). 


e Equation (72) is limited to linearised equations of motion only. However, this is not a problem 
for the spacetimes considered in ref. [45], all of which linearise the Einstein equations, so that 
the linearised double copy can be promoted to an exact statement. 


e The twistor formalism is in principal conformally invariant (at least in those cases in which 
the twistor functions do not involve the infinity twistors for a given spacetime, that break 
conformal invariance). This means that the Weyl double copy should immediately extend 
to conformally flat spacetimes, thus formalising a preliminary observation made in ref. [49]. 
One could, for example, take a given set of scalar, electromagnetic and gravity fields that 
are linked by the Weyl double copy, and conformally transform them directly to a desired 
spacetime, thus achieving a double copy on a curved background. We leave a full investigation 
of this interesting possibility to future work. 


In order for the above to constitute a full derivation of the Weyl double copy, it must be the case 
that all possible vacuum type D spacetimes can be obtained using twistor functions of the form of 
eq. (55). That this is indeed the case has been argued by Haslehurst and Penrose in ref. , and 
general arguments may also be given. Type D vacuum solutions are distinguished by the presence 
of two distinct shear-free null geodesic congruences. All such congruences (in Minkowski space) 
can be obtained as the zero sets of twistor functions, by a result known as the Kerr theorem (see 
also e.g. refs. (125][126}). In the present case, these twistor functions are precisely those appearing 
in the denominator of eq. (55). 


3.2 Example: Schwarzschild & Taub-NUT 


A canonical example is that of the Schwarzschild black hole, which is not (anti-)self-dual by itself. 
However, it is known that duality transformations map out the parameter space of a 
general Taub-NUT solution with Schwarzschild mass M and NUT charge N. Thus, if we restrict 
to the self-dual part of the Weyl tensor only, we will obtain self-dual Taub-NUT with a fixed 
relationship between M and N. To obtain this using the above construction, one may take a family 
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of functions as in eq. (55), where a suitable choice for Qag is 4] 


0 0 0 -1 
Qa =z] oro of (73) 
-1 00 0 
so that the quadratic form in eq. becomes 
y = 717? — 2973 
=w' ty — wry. (74) 


This must be evaluated subject to the incidence relation of eq. (31), where 
, 1 } 
vt ( E aor (75) 
r-twy t-z 
With my given by eq. (56), we then find 


i 


) Va 
TE =7 ET se me E E (77) 


(z+iy) zFr 
and comparing eq. with eq. (57), we find 


(E(x + iy) + 22€ — (x — iy)) . (76) 


This has roots at 


iV2 
(a + ty)” 
From eq. (61), the biadjoint scalar function @ associated with this solution is given by 


N(x) = 


i 
PaT 


This agrees with the function S(x) presented in ref. up to an overall normalisation constant. 
However, their function S(x) is itself only defined up to an overall constant, so this is not a problem. 


(79) 


The field strength and Weyl spinors associated with the electromagnetic and gravity solutions will 
have the form of eqs. (70), where the two principal spinors in the present case are 


t= t= 
‚,=|ą|1 ,= [1 . 80 
av = (22), be = (15) (80) 
We may relate these to the Kerr-Schild double copy of ref. as follows. In a real spacetime, we 


construct the principal null directions corresponding to the real spinors by combining each of the 
latter with their complex conjugates and the appropriate Infeld-van-der-Waerden symbols: 


Ko) — a asod®, KO = Bar Baot“. (81) 


‘We use the self-dual analogue of the anti-self-dual twistor function presented in refs. 142). 
Reference actually presents results for the Kerr solution, but this reduces to the Schwarzschid solution when 
the angular momentum parameter a is taken to zero. 
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Using the spinors of eq. (66), one finds 


kÒ = 5 Ei e a |e); ma 


such that the explicit forms of eq. yield 


ry 2 T Yz 
Ke?) o — (1, a ). (83) 


ra r r r 


The proportional sign here arises from the fact that the spinors of eq. were defined in projective 
space, and thus can be renormalised by a position-dependent factor. This allows us to remove the 
prefactor on the right-hand side of eq. (83), and one recovers the two possible choices of Kerr-Schild 
vector k,, for the Schwarzschild spacetime [44]. 


3.3 Examples of general Petrov type 


Above, we have seen that it is possible to derive the type D Weyl double copy, by choosing ap- 
propriate functions in the twistor space product of eq. (52). However, the spacetime form of the 
(mixed) double copy, eq. (25), is not intrinsically limited to producing type D solutions only. Thus 
the question naturally arises as to whether arbitrary Petrov types are possible. This was briefly 
considered in ref. [120], which presented examples with Petrov types N and III (see also ref. 
for a discussion of Petrov type N). These examples utilised a particularly well-studied class of holo- 
morphic twistor functions, namely elementary states (see e.g. ref. [125]), which consist of ratios of 
factors of the form (AgZ®), where Aa is a constant dual twistor. Elementary states were originally 
intended as alternatives to plane wave states, for the purposes of examining scattering processes 
in twistor space. However, they were recently reconsidered from a different point of view, as the 
twistor functions associated with certain topologically non-trivial electromagnetic fields. Refer- 
ence [143] pointed out that the field associated with the zeroes of the twistor function (AaZ®“) is 
an electromagnetic Hopf knot, any pair of whose electric (or magnetic) field lines are linked. Ref- 
erence [144] generalised this further, by considering the Penrose transform of the family of twistor 
functions 


fr(Z*) = (Aa Z°) (BZP), (84) 


where h is the helicity of the resulting field in spacetime. We may write the dual twistors in eq. 
as 
Aa = (A4, A^), Ba = (Ba, B^). (85) 


Furthermore, we follow ref. in defining the calligraphic quantities 
PaļAaZ°]= A" m4,  pel[BgZ°] =B” Te, (86) 


such that 
A” = ix^ Aa +A”, (87) 


'6In a slight abuse of notation we have used A and B to denote dual twistors, as well as their associated Weyl- 
spinors. However, the nature of the index in each case makes this notation unambiguous. 
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and similarly for B®’. Then the corresponding solutions of the massless free field equations were 


found to be 


2 2h+1 
PAi.. AL, (2) = (ar) Aa, -e AA (88) 


where AAA poe 
, B? AŻ — A*B 
_ B AA’ _ 

Q=ApB”’, yor =i A;B? 
The spin-1 and spin-2 versions of eq. correspond to a null electromagnetic spinor and type 
N Weyl spinor respectively. The general field is referred to as a spin-N Hopfion, given that its 
spacetime topology is related to the well-known Hopf fibration. 


(89) 


A further generalisation was presented in ref. , which considered the family of twistor functions 
(C Z7) =) (Ds Z®)rhre—D) 
(AaZ?)(BpZP horn 


Here h is the helicity as before, and np, ne E€ Z. Defining calligraphic spinorial quantities as in 
eq. (sc), the corresponding spacetime fields via the Penrose transform were found to be 


(AgCl) POD (A py DP )hlne-I) 
(Ap BB’ )h(nptne)+1 


12°) = (90) 


PAs...Ab, (x) = Aa, oe -Aay,)- (91) 


These fields again have a non-trivial topology: in the electromagnetic case, electric field lines 
correspond to torus knots, where nz and np denote the toroidal and poloidal winding numbers. A 
similar geometry (involving gravitoelectric field lines, defined from the parity-even part of the Weyl 
tensor) can be ascertained in the spin-2 case. Again, the electromagnetic field strength is null, and 
the Weyl spinor is type N. Reference sought to generalise this, by constructing gravitational 
solutions with different Petrov types. In particular, the following Penrose transform was noted [7] 


1 N Nab 
(Aa Zo) a (B, Z0) ” [ABH 


Aca Aa Bar ...Ba 


2h) 


Bay, 03 


(ia b b+1 


where 
Nas = (-1)* ( ase i; (93) 


The Weyl double copy properties of torus knots were considered using methods similar to refs. 
in ref. , although the form of the biadjoint field was not explicitly discussed there. How- 
ever, the above results fit very nicely into the twistor space picture for obtaining the Weyl double 
copy. Starting with eqs. [88), one may take the cases with h = 0, h = 1 and h = 2 as the 
scalar, electromagnetic and gravity functions in the twistor space product of eq. (52), such that the 
corresponding spacetime fields are 
a = 2 J Awd ae adod 

o= Ie þa B' = (ar) (AB), PABICID! = (G5) (a'Ap Ac Ap), 


(94) 


17The results of ref. appear not to include an overall combinatorial factor, which we have explicitly instated 
here. 
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from which it is straightforward to verify eq. (25). A similar analysis can be carried out for eqs. 


and (91). 


Next, consider eq. (92). If a = b = 0, one obtains the scalar field of eq. (94). One may construct 
twistor functions of homogeneity —4 by choosing (a, b) = (1,1) or (0,2), leading to the two respective 
electromagnetic spinors 


3 3 
(1,1) 2 (0,2) 2 
/ / = 2p > d / p] / / = Shi he 1 1). 95 
B= -e(a p) AaB a= (Geta) Aaron ae 
Using these in the mixed Weyl double copy of eq. (25), one can generate a number of different Weyl 
spinors: 


2 5 
PRG) = 4 (arr) Ay AB'BoBpn, 


Q|z — 
(1,1)x (0,2) 2 ý 
AED =-2 (ayp) ArArAoBon (96) 
5 
(0,2) (0,2) _ 2 
Py BCID! = (ao) Ay AB Ac Ap), (97) 


which is entirely consistent with the rule for combining the corresponding functions in twistor space. 
The first and third of these examples are Petrov type D and N respectively. However, the second 
(as already noted in ref. (146|) is Petrov type II, thus going beyond the original formulation of the 
Weyl double copy in ref. 


We may go further than the above (and the results of ref. |120)) by seeking Weyl double copy 

examples with Petrov types I and II, where again we may rely on elementary states. However, we 

will see that this leads to a generalisation of the Weyl double copy formula of eq. (51). Let us 

illustrate this with the simpler case of type II solutions. Consider the twistor function 
(11) C27 CF ny 


se (Aaz PBZ (Ata (BT an) 
where Aa, Ba and Ca are constant dual twistors, and we have defined calligraphic spinors as in 
eq. (86). This has homogeneity —6 and thus represents a gravity solution in spacetime, where the 
appropriate Penrose transform can be written as 


(98) 


(1) _ pE (2,3) 
Vypop = C QA B'OID'EN (99) 
where 
(2,3) = 1 TAIT BITC’ T DIT EB! d y 
QA BIC! D! EI ri $ (AC ng) (Bng) Ty 
10 
= [app Adv dc Bo Be, aoo) 


and we have used eq. in the second line. To see why eq. yields a type II Weyl spinor, note 
that one may rewrite the spinorial quantity on the right-hand side of eq. (100) as 


1 
Ay Ap Ac Bp Bry = 5 [BAr Bw + 2Bpr Aqa] Ag'Ac Bpr, (101) 
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such that the Weyl spinor of eq. may be rewritten as 


2 
Uh ee = app 4 AB Be Fon, (102) 
where cA] cB) 
Fal =a } 2 AB) A (103) 


Equation is then manifestly of type II as required, and to show that it may be obtained by 
a Weyl double copy, we must find twistor functions that can be substituted in eq. so as to 
reproduce eq. (98). In fact, we need only use the twistor functions related to the electromagnetic 
Hopfions discussed above. To this end, consider the homogeneity —4 functions related to the two 
electromagnetic spinors (95). These are 


(0,2) _ 1 O 1 
EM (AaZ*)(BgZ?)3  [rA][rB]3 
aD 1 = 1 
EM > ZPB APEEP oe 
where we have used the notations of eqs. [86), as well as the homogeneity —2 function 
1 1 
(0,0) — = : 105 
fO > ZABA FAE] m 
Upon applying the Schouten identity 
[CA] [7B] + [BC][TtA] + [AB] [rC] = 0, (106) 
one may then rewrite eq. as [$] 
1 [CB] [CA] an1 
aD _ (1,1) (0,2) | (1,1) 
grav. — f (0.0) fem ( [AB] fem ' [AB] frm y (107) 
The spacetime field corresponding to eq. (107) is a Weyl tensor of the form 
(I) _ 1 fa [CA] 0,2) aD [CB] a) aD 
Vypgop — $ oh eB = ATAB AB PCD) , (108) 


in agreement with eqs. [103). Here, we have combined the electromagnetic functions to make 
a gravity function in twistor space, and only then carried out the Penrose transform. This correctly 
keeps track of combinatorial factors resulting from the multiplicities of the poles in twistor space, 
and the result of our procedure is that we obtain a generalised double copy formula (eq. (L08)), 
containing a sum of two distinct terms, each with the structure of eq. 


18This factorization is not unique. The other possibility is coe = ol a ( Ga a = Sa ee ), where 


Jer is given in eq. (116). 
19\We have chosen to keep out factors of spinor brackets in eq. (108), but could just as easily have absorbed these 
into the definitions of the electromagnetic spinors on the right-hand side. 
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We may carry out a similar analysis for Petrov type I, by considering e.g. the twistor function 


E’ 2 
pO. = (C* TE) (109) 
grav. (AC na)t(BE' Ty) , 


whose Penrose transform yields 


o aaa cath (110) 
with 
O opam = : - f rE T 
2ri J (AF rer)4 (BH Try)’ 
= -pig AAs Ac Bo Br Ben. (111) 
But 


1 
Ay Ap Ac Bp Be Br = 5 [BA Bry Aca Bp + Ap Ar Biy Bp + Be Bp Ava Ap] Ac Bp, 


(112) 
such that eq. (110) gives 
4 
eee AB [BICA[CB] A Bp + CA}? By Bp + [CBP AAB] Ac'Bpn 
1 
= ——— Ay Bai Dor Ep: 113 
ABE (Bg Do Ep, (113) 
where 
[AC] 2(BC] 
Dor = (3 — V5) — Boa + — Ao, 
[AC] 2(BC] 
Ep = (3 5)—— Bp + — Ac. 114 
D (3+ V5) BD + TAB) c (114) 
Equation (113) is manifestly of type I as required. As for the type II example of eq. (108), it may 
be written as a superposition of pure Weyl double copies. By repeated application of eq. (106), we 
may rewrite the twistor function of eq. (109) as P| 
yl CB] .02) [CA] pa, [C-A] CO 4 [CB] fan (115) 
grav. f (0:0) [AB] EM [AB] EM [AB] EM ' [AB] EM }? 
where we have introduced a third homogeneity —4 function 
(2,0) 1 1 
= = 116 
fM = AZo (BZP) MAPE] ie 
20Similar to the type II case, this factorization is not unique. The other possibility is Moa = 
gam fir (Se (2,0) | [e8]? F002) — gales oie) 
F000) J EM [AB)2/EM ' [ABP IEM [AB]? JEM jJ: 


22 


with the respective electromagnetic spinor 


3 
(2,0) 2 
1p = | s (Bp. 11 


Expanding and transforming each term separately to position space, one finds 


(0) 1 [3 |CAJ[CB] 1 
Vy picrp! = $ Boa aoao — 0 o] 
CBP (02) a1 _ 1 [CA]? 


(1,1) (2,0) 
ABP AB PCD) 2 [AB]? Pagh 688] ` (118) 


As in the previous type II example, this is a sum of pure double copy terms. However, in both 
cases, there is considerable choice in how one presents the final results. Returning to the simpler 
type II example of eq. (108), one may define the alternative electromagnetic spinor 


[CA] (0,2) 4 [CB] 


Pare = 3TABJ AB — TAB] gap oy) 


which is guaranteed to solve the massless free field equation given that the two terms on the right- 
hand side are themselves solutions, and thus may be linearly superposed. Equation (108) then 
becomes 


1 

von = usd, azo) 
which is of pure Weyl double copy form. The reader may be worried that there are apparently dif- 
ferent double copy formulae that can be written down that relate different electromagnetic solutions 
to a given gravity solution. However, this is in fact neither surprising nor profound. The Penrose 
transform used here is limited to the linearised gauge and gravity theories only, as are our examples 
of gauge / gravity solutions, such that the ambiguity in associating a given gravity solution with 
a given pair of electromagnetic solutions is precisely that associated with being able to linearly 
superpose the latter. Notably, the individual terms in eqs. correspond to Weyl spinors 
of restricted Petrov type, such that the superpositions involved correspond to the known property 
that, in the linearised theory, one may superpose solutions to create different Petrov types. We 
have thus succeeded in providing Weyl double copy examples of more general Petrov type, but in 
a rather artificial way. One may therefore question the utility of the twistor approach (and indeed 
the Weyl double copy in general) for these solutions. However, what the twistor framework does is 
provide an interesting way to classify possible double copy formulae, in that the problem of finding 
the different single copies of a given Weyl tensor amounts to obtaining the different factorizations 
of the related twistorial function. It also provides a motivation for why particular solutions may be 
interesting even at linearised level (e.g. the identification of elementary states with Hopfions and 
torus knots [143H146]). It would of course be very interesting to find examples of arbitrary Petrov 
type where exact — or at the very least non-linear — solutions are related. 


3.4 A possible objection 


In the previous sections, we have outlined a derivation of the Weyl double copy, that relies on a 
certain product of holomorphic twistor functions in projective twistor space. However, this should 
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rightly incur the wrath of any sensible twistor theorist: as we discussed in section 2.2} the “func- 
tions” we have discussed above are not actually functions, but representatives of cohomology classes. 
Each spin-n (positive helicity) massless free field in spacetime corresponds to a particular element 
(cohomology class) from the group H! (PT, —n — 2), and the interpretation of the product of eq. (52) 
is then not at all clear 24] 


In more pedestrian terms, the twistor function corresponding to a given spacetime field is not 
unique, but may be redefined by adding functions whose singularities lie on only one side of the 
contour I on the Riemann sphere corresponding to a given spacetime point. Then, the product 
of eq. (52) ) that is needed to obtain the Weyl double copy in position space appears incompatible 
with the ability to perform equivalence relations according to eq. (42), in that the order of these 
operations does not commute. To illustrate this point, it is sufficient to consider redefining the 
twistor functions in the numerator of eq. (52), according to 


ii) = Vea, (121) 


where y(Z) has homogeneity —4, and contains poles either in the northern or southern hemisphere 
when restricted to the Riemann sphere of spacetime point zx, ge ‘ah both. By ea S the 
functions f(Z%) give rise to the same electromagnetic spinors of ye ) (a x) as the functions FO (Z. 
However, forming the product of eq. (52) 52) for the redefined functions leads to the twistor function 


Po EZ). Fo 2 Fe) , Aine. IZ 
iz") Fas) tL Raa) + Aza” ne) 


Both the second and third terms on the right-hand side have homogeneity —6, and thus the right- 
hand side gives rise to a solution of the massless spin 2 free field equation in spacetime. However, 
the second term on the right-hand side involves the original functions FO (Z), and thus will have 
poles in both the northern and southern hemispheres of the Riemann sphere of x. Recognising the 
first term on the right-hand side as our original gravity function in twistor space, we thus see that 
eq. (122) does not correspond to an equivalence relation of the form of eq. (42). Consequently, the 
transformation on the right-hand side will gives rise to a different spacetime gravity solution in 
general. 


If we instead take given representative members of the equivalence class of functions for ( fË, F(Z) 
and form the product of eq. (52), we are indeed free to make redefinitions according to eq. (42). 
That is, the transformations 


1 Q 2 a 1 Q 2 Q 
FZO) IZel) 
f(Z*) f(Z*) 
do indeed yield equivalent gravity solutions. However, we are then faced with the puzzle of how to 


pick out what these representative members are meant to be, given that all possible choices of the 
classes of function entering eq. are meant to be equivalent! 


+ fn(Z°) + fs(Z*) (123) 


21We thank Prof. Edward Witten for comments leading to the present discussion. 


24 


The above puzzle, whilst interesting, does not appear to pose an obstacle to deriving the Weyl 
double copy in spacetime. All one has to do to achieve the latter is to pick suitable representatives 
from each cohomology class, chosen by construction so as to obtain the type D Weyl double copy of 
eq. (51). Put another way, one only needs to verify the following statement: for particular elements 
(cohomology classes) from the groups H'(PT,—2), H'(PT,—4) and H'(PT,—6), a representative 
of each class exists such that the corresponding spacetime fields obey eq. (51). This is a much 
weaker statement than requiring a complete map between the classes themselves i.e. intepreting 
the product of eq. as providing a general map: 


H! (PT, —2) x H'(PT, —4) x H'(PT, —4) > H! (PT, —6), (124) 


which may or may not be achievable. The validity of the weaker statement above is demonstrated 
explicitly in section [3.1] but whether or not anything more general can be said is certainly worth 
investigating, as it is clearly related to central questions regarding the validity and scope of the 
double copy, including to exact solutions of arbitrary Petrov type. 


We note also that from a physics point of view, the situation is highly reminiscent of the well- 
known BCJ double copy for (quantum) scattering amplitudes (1}(2I, in which gravity amplitudes 
are expressed as a sum of terms, each involving a product of kinematic factors {n;} obtained from 
gauge theory amplitudes. These numerators are gauge-dependent, but such that the total amplitude 
is gauge-invariant. The double copy structure is not manifest in arbitrary gauges, and one must 
make generalised gauge transformations (including also field definitions in general) in order to put 
the numerators into a specific “BCJ-dual” form, so that the double copy can be carried out. This 
problem already occurs at tree-level, and if a given set of such numerators is subjected to a gauge 
transformation 


for some 6;, the double copy formula will generate unwanted terms in the gravity amplitude, that 
threaten the gauge-invariance of the latter. It is possible to set up the double copy in a more 
gauge-invariant manner, but at the expense of having to introduce additional correction terms on 
the gravity side, to cancel out the unwanted contributions . Although the situation here is not 
exactly identical (i.e. the equivalence transformations of eq. (42) do not correspond to spacetime 
gauge transformations), it may well be that some similar procedure in twistor space can be defined, 
so that full invariance with respect to equivalence transformations is made manifest. Any such 
procedure presumably faces the additional barrier of having to be interpretable in sheaf cohomo- 
logical terms, but there is again hope. For example, products of twistor space cohomology classes 
have been discussed in earlier literature regarding twistor diagrams for scattering amplitudes (see 
e.g. for reviews). Some of these techniques may be adaptable to the present case of 
classical solutions, and there may also be existing results from the algebraic geometry literature 
regarding maps similar to those required here (although we do not know of anything at the time 
of writing). 


Throughout, we have been discussing twistor cohomology classes using the language of sheaf co- 
homology (or alternatively Cech cohomology, which is an appropriate approximation). However, 
another formulation of the Penrose transform exists, in which the twistor functions become differ- 
ential forms, and are to be interpreted as Dolbeault cohomology classes (see e.g. for a review). 
That this is equivalent to the above approach follows from known isomorphisms between Cech and 
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Dolbeault cohomology groups. It would certainly be interesting to try to reformulate our derivation 
of the Weyl double copy in the Dolbeault approach, as this is clearly related to whether the double 
copy has a genuinely twistorial interpretation. 


3.5 The Weyl double copy for anti-self-dual fields 


In the previous sections, as in ref. [120], we have addressed the Weyl double copy for self-dual fields 
i.e. those with primed spinor indices. In this section, we extend this discussion to anti-self-dual 
fields. As reviewed in section [2.2] there are two Penrose transforms one may consider for anti-self- 
dual fields. The first (eq. (43)) simply consists of replacing twistors with dual twistors, and it is 
straightforward to see that the derivation of the type D Weyl double copy in terms of anti-self-dual 
fields proceeds similarly to the case of self-dual fields discussed above. That is, one may consider 
the family of functions 


j= = [Q wawa] E (126) 


with Q® a constant matrix. In the Penrose transform, this is to be evaluated subject to the 
incidence relation of eq. (36), and one may choose homogeneous coordinates 


Aa = (1,7) (127) 


such that the quadratic form appearing on the right-hand side of eq. (126) may be written as 
X= pr [QÊWa We] = Ñ= E)n — m(2))(n — mo(2)), (128) 


for some spacetime-dependent functions N and 7. Carrying out the Penrose transforms for m € 
{1,2,3} yields spacetime fields 


N(a Ñ? (ax NB8 (a 
an M ( L, QAB = -re aay ABCD = rE a(abnacbp), (129) 


obeying the Weyl double copy formula of eq. (23). 


S 


One may also consider using the (non-dual) twistor space Penrose transform of eq. (44), but the 
complication then arises of how to form a product in twistor space (i.e. before or after the derivatives 
are applied). In section|3.1| each quantity entering the twistor space product must be interpretable 
by itself as corresponding to a spacetime field, after restriction to a given spacetime point. In 
eq. (44), the restriction to a given spacetime point happens after the function f(Z°) has already 
been differentiated, which suggests that we define a twistor-space product in terms of differentiated 
quantities: 


o oð o 


(Z°). (130) 
A twistorial double copy for anti-self-dual fields can then be written as 


FABÍCD) 
ABCD 7 F ` (131) 
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We can at least show that such a relationship holds in particular cases. For example, a suitable 
function to be entered into eq. for the (anti-self-dual) Coulomb solution is [140] 


J™ = log ($) , (132) 


with 
(133) 


From eq. (29), we then find 


fô) 
apa and Awad zQ =0, (134) 
so that it is straightforward to compute 
eM _ QAQB _ Od 
SAB el Q? , QA = 3wvA Q. (135) 


The anti-self-dual Schwarzschild / Taub-NUT solution can be obtained from the following twistor 


function for use in eq. (141): 
Q 


je = 190g (136) 
with Q and P as defined in (133). Then 
perav QaQBQcQpd 
ABCD = ~~ Qs (137) 


Finally, comparing equations (135) and (137), we see that the double copy formula (131) is indeed 
verified, with 


Jeg (138) 
We remark that the expression for the scalar twistor function is the same as that used for the self- 


dual analysis of section [3.2] as must be the case. Furthermore, different choices of the quadratic 
form Q (subject to the conditions of eq. (134)) will map out the space of type D vacuum solu- 


tions [139]. 


It is possible to extend the above to general families of solutions. Firstly, recalling the definition of 
Z® eqn. (30) 
Ze = (Z°, Zt, Z2, Z8’) = (w®, wa, ny) = (w4, Ta) (139) 


we notice that 
Q = Qag Z“ ZÊ = Z! 2? — Z0 Z? (140) 


was chosen exactly such that 
o 
Qa = ə zQ = (T1, To) = (em) a, (141) 
w 


so that Q4 is just my rotated by a Levi-Civita symbol. It is then straightforward to show that we 
can write 


P:[Q] = (AW ry) (B? mB) (142) 
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with A^ and B®’ defined as in (87). Then, using (141) and (142) we can write the integrand for 
the Coulomb solution of (132) as: 


ð ô seml l 
aa | = — (em) (em) B (AT ry) (BB Tg)? (143) 
= — (en)a(em) a fii? 


with poy the function oa (SO) in the self-dual transform (see eq. (104)). Similarly, for the 
Schwarzschild solution of eq. we have 


1 
(A^ Ta)’ (BP Tp)’ (144) 


=2(er) (em) B(em)o(em) rfp? 


a 


Bw4 Aw Aw! AwP jem =2(er) a(et) B(ew)c(eT™)p 


Above we described the double copy for Type D. In order to progress to more general families, we 
will first use the results above to find the f©™’s which map to f ea and Har defined in (104) and 
(116). Making the ansatz 


~(0,2) Q Q iā 
EM = Fale os (2) (145) 
with Q and P as before and 
R= RaZ% = R” ry, (146) 
we have a9 ə QQ 
7(0,2) _ QAQB 
DWA JB EM = RQ En 
Then 
pe | sagen iii | = eater) . 
DAD R = A B 7 7 7 
Jwi Ju EM (AT ry )(BB'Tg) (RETA)? (148) 
A! _ gA 
AE (emalem)a fiir 
Similarly, we have 
-(2,0) Q Q 
ee ge os ($ ) (149) 
with Q and P as before and 
S = Sa Z% = Sry, SX =A”. (150) 
Finally, the anti-self-dual analogue of the gravity function (109) will be 
2 2 3 
z(a Q Q Q Q 
Faw. = [AQ +92% + Gz + aa log ($) (151) 
with P, Q, R, S defined as before. If we choose 
CAJ][CB cB]? CA]? CB]|[CA 
Gi =-14 oho, G2 = +B Gs = 3 Gar. ga = Se (152) 
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then the double copy factorisation proceeds by direct analogy to (115) and the subsequent discus- 
sion. 


The double copy formula of eq. is perhaps less desirable than the form based on dual twistors, 
in that it ceases to be a simple product, and thus appears to offer no additional advantages with 
respect to the spacetime double copy formalism. Note also that the same objections regarding 
how to interpret the procedure in cohomological terms apply here. The twistor “functions” to 
be entered into eq. are actually cohomology classes, which in this case are elements of the 
sheaf cohomology group H!(PT, O(n — 2)), for a spin n field. Differentiating 2n times maps each 
cohomology class into an element of H!(PT, O(—n — 2)), similar to the case of self-dual fields. Once 
again, we may take the pragmatic view that in order to generate a particular spacetime double 
copy, it is sufficient to show that particular representatives of the cohomology classes may be found 
in twistor space, that achieve the desired spacetime relationship. 


4 Conclusion 


In this paper, we have examined the Weyl double copy that relates solutions of biadjoint scalar, 
gauge and gravity theories, using a twistor-space formalism initiated in ref. [120]. The latter argues 
that each instance of the Weyl double copy in spacetime can be associated with a certain product 
of functions in twistor space. We have provided full details of how this formalism is sufficient to 
derive the previously noted form and scope of the Weyl double copy, namely the fact that it applies 
to arbitrary vacuum type D solutions. We have also gone further than ref. in providing exam- 
ples of Petrov type I and II solutions in gravity, in addition to types III, D and N. However, such 
solutions are limited to linearised level, which is ultimately due to the limitations of the Penrose 
transform itself. We have also shown how similar arguments can be used to derive spacetime double 
copy formulae for anti-self-dual fields, as well as self-dual ones. 


Care must be taken in how to interpret the twistor space double copy, given that it apparently 
involves multiplying together twistor functions. In the Penrose transform, the “functions” are in 
fact cohomology classes (i.e. elements of sheaf cohomology groups). Deriving a given instance of 
the Weyl double copy then amounts to showing the existence of appropriate representations of each 
cohomology class, such that the functions entering a particular instance of the spacetime Weyl 
double copy are indeed related by a twistor-space product. This is a far cry from demanding a 
map between the relevant cohomology groups themselves, and the investigation of whether a more 
rigorous twistor-space interpretation exists deserves further investigation, as it may shed further 
light on the ultimate origins and scope of the double copy itself. It may also open up the possibility 
to look at fully non-linear solutions. Work on these issues is in progress. 
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